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INFLUENCE OF SHEAR DEFORMATION OF THE CROSS SECTION ON 
TORSIONAL FREQUENCIES OF BOX BEAMS 
By Edwin T. Kruszewskl and William W. Davenport 

SUMAEI 


An exact analysis has heen carried out on the torsional vibrations 
of a four— flange box beam with cross sections which can change shape 
beca^lse tbe stiffness of the bulkheads is finite. The effect of shear 
deformation of the cross section on the torsional frequencies is illus- 
trated by numerical calculations . An approximate method for quickly 
estimating" the effects of bulkhead shear stiffness on the torsional 
frequencies of box beams has been devised. 


INTRODUCTION 


In an experimental investigation described in reference 1, analyses 
based on the assunptlon that the changes in the shape of the cross sec- 
tion are negligible were foimd to be completely inadequate in predicting 
the experimental torsional frequencies of a thin-walled tiibe of rectan- 
gular cross section. One form of cross-sectional distortion - due to 
local deflections normal to the surface of the covers and webs - was inves 
tlgated in reference 2. The magnitude of this effect on torsional 
.frequencies, however, was shown to be practically negligible for the 
particular test beam and frequency range investigated in reference 1. 
Another type of cross-sectional distortion - due to overall changes of 
shape of a sheeiring nature - coiold be an important Influence if the bulk- 
heads were not rigid. Although the effect of bulkhead shear deformation 
has been considered in some analyses (see, for example, ref. 5 ) ^ no 
investigations of the importance of this effect on torsional frequencies 
seem to exist. 

0?he importance of the shear stiffness of bulkheads in torsional 
vibrations is assessed in the present paper by means of an analysis of 
a four -flange box beam containing bulkheads with finite shear stiffness. 
Frequency equations are derived for torsional vibrations of a uniform 
free-free beam, and numerical results obtained by use of these equations 
are shown. A set of ctirves ’based on an approximate solution, from which 
the effects of cross-sectional shear deformation can be quickly estimated, 
is also presented. 
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SYMBOLS 


A, B 

C 

E 

G 

Ge 

J 

K 

L 


area of flange 

parameters defined in equations (22) 

inertial coupling constant. 

Young's modulus of elastici’ty 
shear modulus of elasticity 
effective shear modulus of bulkheads 

mass moment of inertia per unit length about y-axis 
mass moment of inertia per imit length about z-axis 
mass polar moment of inertia per unit length, ly + 

torsional-stiffness constant 

restralnt-of -warping parameter (see eqs. (22)) 

half-length of beam 



M parameter, — ^1 - 

P freq.uency parameter for special case (see eq. (A9)) 

S bulkhead stiffness parameter (see eqs. (22)) 

T maximum kinetic energy 

U maximum strain energy 

aQ,a^^,b^,c^ Fourier series coefficients 

a half -depth of beam 

b half -width of beam 
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i, m, n 



X, y, z 


integers 

frequency coefficient, 
frequency coefficient for uniform shear mode 
thickness of cover sheet 

thickness of weh 

displacement of flange in x-dlrectlon, y-direction, and 
z -direction 

coordinates defined in figure 1 
longitudinal strain in flange 



’'C 

’'h 


shear strain in cover 
shear strain in web 
shear strain in bulkhead 


e 


X 



average angle of twist 

natijral torsional frequency of four-flange box beam 
frequency of cross-sectional imlform shear mode 

mass density of the actual beam 
Kronecker delta (lifi=Oj 0 Ifl^O) 


Subscripts i 

1, n integers 

m mode number 


rig 


rigid bulkheads 



k 
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THEOREOICAL ANALYSIS 
Basic Equations 


In order to investigate the influence of shear stiffness of bulk- 
heads on the torsional frequencies of- uniform box beams, it is conven- 
ient to make an idealization of the actual structure. In reference 4, 
the idealization of the webs and covers of a beam into flanges and a 
shear -carrying sheet was successfully used in an analysis of a beam with 
rigid bulkheads. In the present analysis, the same idealization is used 
for a beam which contains nonrigid bulkheads. For simplicity of calcu- 
lation, the cross section of the four-flange box beam is assumed to be 
doubly synmetrlcal. 

Assumptions .- In the present analysis the following assunptions are 

made: 


(1) The flanges of the beam carry only normal stresses. 

( 2 ) The sheets connecting the flanges carry only shear. 

( 3 ) The beam contains continuously distributed. Independently acting 
bTilkheads that have finite shear stiffness in their planes but are 
entirely free to warp out of their planes. 

(4) The influence of longitudinal inertia is negligible. 


In accordance with the foregoing assunptions and the double symmetry 
of the cross section, the displacement of any point on the cross section 
can be defined in terms of u, v, and w, the displacements of one of 
the flanges in the x-, y-, and z-directions, respectively. (See 

fig. 1.) 


Strain relations .- The longitudinal strain in the flanges e^, and 
the shear strains in the covers, webs, and bulkheads 7y, and 

respectively) are given in terms of the displacements u, v, and w 
as 


^ dx 


( 1 ) 


y = ^ + U 

C dx b 


(2) 
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„ _ dw u 

^ dx a 


( 3 ) 




V _ w 

a "b 




where a and h are the half -depth and the half -width of the beam, 
respectively. 

If the bulkheads were assumed to be rigid in their own planes, v 
arifi w would each be proportional to the twist of the cross section in 
such a manner that 7 ^ = 0. Ihe distortions of the structure could then 
be defined as in reference in terms of u and the twist of the cross 
section. Since the bulkheads have a finite shear stiffness, the cross 
section is allowed to distort and v and w are considered separately. 
It is then convenient to define an average twist 6 as 


e 



(5) 


From equations (!{•) and (5) , the displacements v and w can now be 
written in terms of 0 and as 



Energy rf>l n.tions . - For the fo\jr -flange beam vibrating in a natural 
mode, the maximum strain energy U is 



where L is the half-length of the beam, t^ is the cover -sheet thick- 
ness, t^ is the web thickness, Aji is the area of a flange, and Gg 
is the effective shear modiolus of the bulkheads . 


The first two tenns of equation (8) represent the energy due to 
the shear strains in the cavers and webs of the beamj these are the only 
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kinds of energy considered in elementary theories. The third term 
represents the contribution of the restraint of warping. The last term 
represents the contrihutioai of the shear strain in the bulkheads . 


The strata energy expressed by equation ( 8 ) can now be expressed 
in terms of u, 7 ^, and 0 as 

" = 20bt, /_^(a g . I gb , ^ - 1 ^ ^ 


dx + 2 Ggab 


r^^dx 


(9) 


In the calculation of the maximum kinetic energy, the inertial 
prqperties of the actual beam are used and the following assumption is 
made: At any point in the cross section, the displacement v in the 

y-directlon is proportional to z, and the displacement w in the 
z-direction is prcportional to y. (See fig. 1.) Thus the maxlmxmi 
kinetic energy can be written as 




dz dy dx 


( 10 ) 


where 03^ is a natxrral torsional frequency, p is the mass density of 

any point in the actual beam, and R is the region in the plane of the 
cross section containing all of the mass elements of the beam. 


When equations ( 6 ) and ( 7 ) are substituted into eqmtlon (lO) , the 
expression for the maximum kinetic energy can be written in terms of 7 ^ 
and 6 as follows: 



(11) 


where 


C 


Iz - ly 


(12) 


The quantity Ip is the mass polar moment of inertia per unit length, 
and 3^ and Ig are the mass moments of inertia per unit length about 
the y- and z-axes; these quantities are, of course, based on the 
assumption of unifoim spanwise mass distribution in the actiial structure. 
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T 


The constant C can be looked upon as an inertial coupling coefficient 
between the rotational and cross-sectional shear motions. 


Method of Analysis 

The following analysis is based on the principle that a natural 
mode of vibration must satisfy the variational equation 

8(U-T) = 0 ( 15 ) 

where the variation is taken with respect to the distortions u, 0, 
and Tjj eind where these displacements satisfy the geometric boundary 

conditions. Application of the techniques of the calculus of variations 
to equation (l 3 )j with Uj 6 , and 7 .^^ as the independent variables, 

yields the differential equations and natural boundary conditions given 
in the appendix. It is more convenient, however, to carry out the exact 
solution of equation ( 13 ) by means of the Rayleigh-Rltz procedure, 

Syrtmetrlcai vibrations of a free-free beam .- Appropriate assuniptions 
for the distortions u, 0, and of a free-free beam in a synmetrical 

mode of vibration are 

CO 

u = sqx + ^ sin (llj.) 

n=l,2,5 


00 

0 = ^ cos 3^ ( 15 ) 

CO 

7-b = X Cji cos ( 16 ) 

n= 0 ,l ,2 


The choice of the particular trigonometric functions used in equa- 
tions (l4) , ( 15 ) , and ( 16 ) was guided by consideration of the orthogonality 
required for the simplification of the strain-energy expression. The 
linear term a^x is Included in the expression for u in order to allow 
the deflection of the tip of the beam to be unrestricted. 
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Substitution of equations (l 4 ), (15), and (16) into equations (9) 
and (11) yields, after subtraction, 

U - T = asbto / ^ ^ ^ Cntelsin^t 

^-l 1_ n=0,l,2 ^ -2 n=^,2 ^ 

00 ~ 2 

r ^ slq sin £S2£ + ^ x dx + 

° n=l,2,3 ^ ^ _ 

2Gat^ r -b X b^Msin^+J X c„fM\sin M2£ _ 

I dx+2EApr^ao + 

n=l,2,3 ^ ^ J -L _ 


J ^f^)cos 2 p dx + 2 G ab T f X cos ^ 

n=l,2,5 ® -'-LVn=^,2 “ 

/ L/ 00 \2 

( X ^ cos ^ j dx - 
L\n= 0 ,l ,2 ^ j 




I W'^/ 


X ^)( X 

•LVn=0,l,2 ^ /\n=0,l,2 


c COB dx 

n T. , 


Differentiating U - T with respect to a^, the a^'s, the b^^'s, 
and the c^'s and setting the respective results equal to zero yields 

n=§,5> ^ ^ ^ ^ ■ 2 + K^)aQ = 0 (I 8 ) 
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^ 

k ^ + A(i«)bi - 1 B(l^)ci - ^ B ^^0 = ° 




(1 = 1,2,3, • • •) 

(19) 

^ A(i*) ^ - 
a D 

+ 8gj^) - B{l>t)^ 

h + [jc%^(l + 60 i) - 


|A(ijt)^Ci . 

. ^ A(1 - 5oi) (.-1)\ 

= 0 (i = • • •) 

(20) 

1 L B(lrt) ^ 

2 a h 

- [i “^^(1 + soi) - 1 

A(ijt)^bi + + &oi) - 


B(iit)^ - S^l 

^oi)]=i - i 

" ^0i)^0 “ ° 


(i = 0,1,2, . . .) 

(21) 


where 
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The imlm nw-n natural frequency is contained in the parameter 
which is in common use in torsional vibration analyses. The quantity J, 
which appears in the ejjpression for ijji, is the well-known torsional 
stiffness constant. The parameter K is associated with the effects 
of restraint of warping, whereas the parameter S is associated with 
the effects of bulkhead stiffness. The parameters A and B are 
siniply geometrical properties of the four -flange box beam. 

Examination of equations (l8) to ( 2 l) shows that the coefficients bQ 

and Cq appear only in equations (20) and (2l), and then only for i = 0. 
When 1 = 0 , equations ( 20 ) and ( 2 l) reduce to 


2kqPhQ - Ck^^CQ = 0 


2CkQ,%Q - - s)cq = 0 


( 23 ) 


The condition for a nontrivial solution of bQ and Cq gives the 


frequency equation 



( 24 ) 


From equations (15) and (16), it can be seen that a given value of 
bQ corresponds to a rigid-body rotation of the beam, whereas a given 
value of Cq corresponds to a bvQihead shear distortion which is uniform 
along the length of the beam. Equation ( 24 ) , therefore, gives the fre- 
quency coefficient for a uniform cross-sectional shear mode. 



(25) 


in addition to the frequency coefficient kq, = 0 for the rigid-body 
torsional mode. The frequency coefficient for the uniform cross- 

sectional shear mode is shown to be a pertinent parameter in the evalua- 
tion of the Inflvtence of the shear stiffness of bulkheads on torsional 
frequencies . 

Wow consider the remaining eqiiations (18) to (2I) . Substitution of 
equation (19) into equation (18) gives 
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ao + ^ 


n= 



(-l)“(njt)a^ = 0 (k2 / 0) ( 26 ) 


By solving equations (l 9 ) j ( 20 ) , and ( 2 l) simultaneously for 
in terms of sq and then suhstituting the results into equation (26), 
the following equation results if it is required that Sq is not equal 
to zero: 


1 + 2 




n-1,2,3 + 




(«}2 . 

(¥? i 0 ) ( 27 ) 


where 



(28) 


Equation (27) is the frequency equation for the symmetrical torsional 
vibratlon of a four-fleinge hox beam, when the Influence of bulkhead shear 
deformation is included. 

Antisymmetrical vibrations of a free-free beam .- For a free-free 
beam vibrating in an antisymmetrical torsional mode, appropriate assump- 
tions for the distortions are 


u = a^ + Sq cos ^ (29) 

n=l,3,5 ^ 

e = J <50) 

n=l,5,5 ^ 

^ c^sin^ (51) 

n=l,5,5 

As in the case of the synmetrical modes of vibration, the choice of 
the particular trigonometric functions was guided by considerations of 
the orthogonality required for the singjllflcatlon of the strain-energy 
expression. The constant term aQ in equation (29) is necessary to allow 
sufficient freedom for the beam to warp. 

By substituting the expressions for u, 9 , and y-^ (eqs. (29) 
to (51)) into equations (9) and (ll) and by following a procedure similar 
to that described for the symmetrical modes of vibration, the following 
frequency equation for the antisymmetrical mode of vibration is derived: 
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t. 








g )St^ ^^Tspi^ l 


n- 1 , 5,3 


1 - C2 






(kS/0) 


(32) 


DiscTission of Limiting Cases 

Before proceeding with a ntimerical evalmtion of the influence of 
the shear stiffness of hulMieads on the torsional vibrations of hox beams, 
a discussion of some limiting cases of the freq^uency equations is desirable. 

The Influence of cross-sectional shear distortions on torsional fre- 
quencies depends on the frequency coefficient of the uniform cross- 

sectional shear mode. The stiffen the cross sections of a beam, the 
higher the frequency of this uniform shear mode. When the cross sections 
become rigid (that is, when Gq enidL, consequently, are infinite), 

the frequency eauatlons can be put into closed forms Identical to those 
in reference k (where the infl\aence of cross-sectional distortions is 
considered negligible) . 

The parameter K is associated with the influence of restraint of 
warping and, if this effect is to be neglected, it is sufficient to set 
K equal to zero. The frequency equations (eqs. (27) and ( 32 )), however, 
become indeterminate for K = 0 . A solution for this special case is 
presented in the appendix and has a particular in?>ortance which will be 
discussed later in the section entitled "tfethod for Estimating Effects." 

Another special limiting case which is of Interest is the box beam 
for which both C and A are equal to zero as, for example, a square 
beam with eqvial wall and cover thickness. For this case, equation ( 20 ), 
for the condition b^ f 0, yields the elementary torsion^ frequency 
eq-uatlon. whereas solution of the remaining equations (eqs. (18), (l9)# 
and (21) ) , which no longer contain b^, results in a frequency equation 
for cross-sectional shear modes (that is, modes which contain only 
and u displacements) . The frequency equations for these shear modes can 
be obtained from equations (27) and (32) by setting both A and C equal 
to zero. 


NUMERICAL RESULTS 


In order to evaluate the influence of shear stiffness of bulkheads 
on torsional vibrations, the frequencies of the rectangular tube shown 
in figure 2 were calculated frcsn equations (27) and (32) for various 
values of kj^ . The tube was assumed to have a width-depth ratio b/a 
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of 5.6, a plan- form aspect ratio L/b of 13 - 3 » ®- thickness ratio 

t;^/tc of 1 . 0 . A value for E/G of 2.65 was also assumed. (These pro- 
portions correspond to those of the tube for which experimental results 
were reported in ref. 1 .) The cross-sectional area of the fla n ges of the 
four-flange box beam was taken as equal to one-sixth of the sectional 
area of the walls adjacent to each comer. The cover and web thicknesses 
of the four-flange box beam were set equal to the wall thickness of the 
tube. The results of these calculations for the first five free-free 
modes (three antisymmetrical and two symmetrical) are shown in figures 2 
and 3. 


In figure 2 , the frequency coefficient k^^ of the four -flange box 

beam is plotted as a function of the frequency coefficient kjj^ for the 

uniform cross-sectional shear mode. Two sets of curves are shown; the 
solid c\irves represent the frequency coefficients obtained from equa- 
tions (27) and (32) when the influence of cross-sectional distortions is 
included, and the dashed curves represent the frequency coefficients 
obtained from a solution of the foxir -flange box beam when the cross- 
sectional distortions are assumed to be negligible. The solid curves 
representing the first five free-free modes are given for values of kjj 
from 2.73 to 25. ° 


The differences between the solid and dashed curves in figure 2 , 
show that the reductions in torsional frequencies due to bulkhead shear 
flexibilities can be of considerable Importance, especially for small 
values of ^q* a tube with fcjjQ = 3«08, for example, the reduction 

in torsional frequencies due to cross-sectional flexibilities, as 
obtained from figure 2 and verified experimental ,ly in reference 1, ranges 
from 16 percent for the first mode to 63 percent for the fifth mode. 

This value of kt>Q Is for the test beam used in reference 1 . The beam 

contained no discrete bulkheads of any kind and therefore depended on 
the bent action of its own walls for cross-sectional stiffness. If, 
however, weightless bulkheads of the same thickness and shear modulus 
as the walls of the beam were spaced at a distance of twice the chord, 
the average value of k^Q would be approximately 22.0, and the reduc- 
tions in torsional frequencies would range from 1 percent for the first 
mode to 8 percent for the fifth mode. 


In figvnre 3 ^ the results of figure 
frequency ratios ^ i 

the frequency of the mth nxde of the four -flange box beam with rigid 
cross sections. It should be noted that the ratios of the frequencies 
and oijj //o>Ti ^ . are identical to the ratios of the 


2 are replotted in terms of the 


(%): 


rig 


is 


‘%./(“^)rlg 

c orresponding 


’o/(^*^)rig 

frequency coefficients 


kgj 


ny^(^m)rig 


and kjj 


b/(^^^0rig 
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0 ?he curves fcxr each, mode in figxare 3 fall so close together that it 
should he possible to draw one curve which would be representative of 
all modes considered. It would seem probable, therefore, that, with 
some simplification of the analysis, a relationship between cd^ i 

and ^ could be obtained that does not depend on the mode. 

Such a relationship would be useful for estimating the effects of cross- 
sectional shearing. 


MEOIBDD FOR ESTIMATING EFPECOB 


It is shown in the appendix that, when the restraint of warping is 
neglected (that is, when EAj. and, ccaaseq.uently, K are set eq.ua! to 

zero), the following frequency eq.uation will result: 



(55) 


Equation ( 33 ) is similar to the equations obtained in reference 2 in 
which the influence of cotq)ling between overall torsion and certain 
cross-sectional or panel vibrations was investigated. The particular 
type of cross-sectional vibration considered in reference 2 , however, 
is one in which the comers of the cross sections do not move with 
respect to each other; thtis, no cross-sectional shear distortions are 
allowed. __ 

Equation (35) is a quadratic equation in | 1 and will 




rig 


yield two real and positive values of S . Only the smaller of 

(%)rlg 

the two roots, however, is of interest. Results of this approximate 
solution are shown by the dashed curve in figure 3 for the particular 
beam considered in figures 2 and 3 and are in good agreement with the 


“^0 


results of the exact solution, even for the lower values of 
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This frequency equation (eq.. (33) )> therefore, can he used to 
estimate quickly the effects of cross-sectional shear stiffness on the 
torsional frequencies of box beams. Solution of the equation for a range 
of values of JyjT^ (which determines the value of C) from 1.0 to 0.05 

is shown in figure i|-. A value for ly/lz of 1-0 corresponds to a 

uniform-walled tube with a width-depth ratio b/a of 1.0, whereas a 
value for ly/lz of O .05 corresponds to a tube with b/a approxi- 
mately equal to 7«0» curves in figure k show that the effect of 

cross-sectional flexibility increases not only with a decrease in 
uniform cross-sectional shear frequency but also with an Increase in 
width-depth ratio. 


Before these curves can be used, the ratio ly/lg 


and the frequency ratios ^'bo/(^): 


most be known 
must be 


evaluated. A value for the frequency coefficient 


V -Tnyrig 


may be 


obtained from any torsional analysis in which the cross sections are 
assumed to be rigid. A value for the frequency coefficient kj^ may 


be obtained from equation ( 25 ), provided an appropriate value of Ge 
for a representative beam cross section is known. In determining Gg 
for an actual structure where the bulkheads contribute most of the cross- 
sectional shear stiffness, the defining relation is 


°e 

G 2LG 




where N is the number of bulkheads and 


(Gtb). 


is the effective shear 


stiffness of a bulkhead. If the bulkheads are solid ^heets, the shear 
modulus of the bulkhead material should be used for G and the bulkhead 
thickness shoiild be used for t^,. For other forms of bulkheads, estimates 
of the effective shear stiffness nnist be made. For beams that contain no 
biilkheads, special consideration, such as that given in appendix B of 
reference 1, is necessary. 


CONCL-UDING RBMAKKS 


The influence of cross-sectional or bulkhead shear deformations on 
the torsional frequencies of box beams has been obtained by means of an 
analysis of a four -flange box beam. For conventional constructions where 
bulkheads are spaced approximately 1 chord apart, the influence of cross- 
sectional shear deformations generally is negligible. For beams that 
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depend on the bent action of their own walls for most of their cross- 
sectional stiffness, the effect of cross-sectional shear deformations is 
considerable. Since the trend in wing design is toward structures with 
fewer bulkheads and higher width-depth ratios, the Inf luence of cross- 
sectional flexibilities on torsional fTeqxiencies could become important. 
Curves based on an approximate solution ar6 presented to allow q.vilck 
approximation of the reduction in torsional frequencies of box beams 
having flexible rather than rigid bulkheads. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 27 > 1955 * 
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* 

APPENDIX 


* 

APPROXIMA.TE SOUJIION 



Al+.^migh a solution for K = 0 can be obtained directly from equa- 
tions (l 8 ) to ( 21 ), a solution based on tbe differential equations mi^t 
be of more interest. 


From eqTiations (9), (ll), and (13), where the variation is taJsen 
independently with respect to u, 0 , and the following differen- 

tial equations and the natural boundary conditions associated with the 
vibration problem may be obtained: 

m 

d ^0 1 , A du ^ 0 . 1 

= - 2 * ^2 - at to ^ 1,2 ® 2 1,2 ’V. 

(Al) 

• 

J^^abto- 1,2 -S)^rb -0 

(A2) 


- - A — + i B + — u = 0 

ab dx 2 clx ab 

(A5) 


k 2 ^ 5u ^ =0 
^ -L 

(Ai^) 


/'sSa-lAi^-^-.^se'' .0 

y dx 2 dx ab j 

(A5) 

* 

(a — - - B \ 1 = ° 

dx 2 dx ^ J ^\-Jj 

(a 6 ) 

% 

The parameters used in these equations ane defined in equations 

( 22 ). 
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In order to neglect the influence of restraint due to warping in 
equations (a 1) to (a 6), it is sufficient to set equal to 0. Ihese 
equations may then be reduced to the following differential equation and 
boundary conditions: 



^ + P^e = 0 
dx^ 

(AT) 


L 

M 60 =0 

^ -L 

CD 

where 

l-02 

(A9) 

For a free -free beam, 

, equation (a 8) will yield 



M ^ - 0 

dx _L ax L 

(AlO) 

The general solution 

for equation (AT) is given by 



0 = sin Px + D 2 cos Px 

(All) 


where and D2 ai^ constants of integration. 

Now, by use of equations (AlO) and (All) and the condition for a 
nontrivial solution for D]^ and D 2 , the following frequency equation 

for a free-free beam may be obtained: 


sin PL cos PL = 0 


(A12) 
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Equation (A12) is satisfied by the relation: 


p2l 2 = (m = 0, 1, 2, . . . ) (AI5) 


By substitution of the definition of p2 from equation (A9) > equa- 
tion (AI 5 ) niay be written as 



S 


kQi' 


1 - 


- l=br 



(m = 0, 1, 2, . . . ) (A14) 


For m = 0, eqmtion (A14) yields the frequencies of the rigid-body 
torsion mode and the uniform cross-sectional shear mode (see eq. ( 25 ))* 
If the case in vhich m = 0 is neglected, equations (A14) may be 
written as 


(^m)rig 



P 


2 


,2 

1 


(“^m)rlg 


^°^m)rlg 



K) 


rig 


1 


2 


1 - C 

(“>Tm)rig 


(“Tm)rig 


= (= 




(m — 1, 2, 5# • • •) 


(AI 5 ) 


At this point, it should be noted that reduces to the 

mth elementary frequency coefficient for the case in which = Oj 

that is, (l%i)rig " where m = 1 , 2 , 5 , . . . 

Thvis, equation (AI 5 ) may be written as 



4 

~ ^0 

2 

1 



2 

“bO' 

(°^m)rig 


("^)rlg 

^ 1 - 

J 


(“Tm)rig 

T 

(“^m)rig 


(AI 6 ) 
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